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Abstract 

Neutrino telescopes with large detection volumes can demonstrate whether the current indica- 
tions of neutrino oscillation are correct or if a better description can be achieved with non-standard 
alternatives. Observations of contained muons produced by atmospheric neutrinos can better con- 
strain the allowed region for oscillations or determine the relevant parameters of non-standard 
models. We analyze the possibility of neutrino telescopes measuring atmospheric neutrino oscilla- 
tions. We suggest adjustments to improve this potential. An addition of four densely-instrumented 
strings to the AMANDA II detector makes oscillation observations feasible. Such a configuration 
is competitive with current and proposed experiments. 
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I. INTRODUCTION 



Experimental observations (Super-Kamiokande |], Kamiokande 0, 1MB [Q] and Soudan 
0]) for atmospheric muon and electron neutrinos have found that the ratio of the number of 
these neutrino species does not agree with theoretical prediction. All experiments measuring 
the flux of solar neutrinos observe a deficit compared with the solar model predictions ||. 
The ratio of muon to electron events observed in atmospheric neutrino interactions is mea- 
sured by most experiments to be less than expected from models of cosmic ray interactions 
in the atmosphere [El. Neutrino oscillation can explain these results. The measurement of 
the up/down asymmetry of this ratio by the Super-Kamionkande collaboration [[]]] is gener- 
ally considered to be the strongest evidence for neutrino oscillations. As oscillations would 
most likely imply that neutrinos have mass, many researchers have fitted the available data 
for the mass difference among different neutrino species. 

While the cumulative evidence for neutrino oscillations is very striking, a definitive proof 
that the observed anomalies are actually due to neutrino oscillations is still missing. In par- 
ticular the current observations of atmospheric neutrinos are consistent with the hypothesis 
of maximal oscillations, but do not exclude some alternative unconventional explana- 
tions, such as neutrino decay [0, micron-scale extra dimensions || or quantum decoherence 
in propagation ||. 

Cosmic rays traveling to Earth will interact with its atmospheric nuclei. These interac- 
tions produce hadrons which decay eventually into atmospheric neutrinos. Up to neutrino 
energies of about 100 GeV, the main result from these interactions is pion production 7r + (7r~). 
These will decay into /i + + + V^) almost 100% of the time. The secondary 

will decay into e + + v e + z/ M (e~ + V e + u^) which will give a + to v e + V e ratio (r) of 
2:1. 

However, depending on the cosmic ray energy and where it interacts in the Earth atmo- 
sphere, the muon might not decay before reaching the ground and corrections to this ratio 
have to be included. Another effect to be included is the Earth's magnetic field action on 
lower energy particles. These effects have been modeled and simulated and one looks for a 
ratio (R) that relates r exp from experimental data to r sim from simulations (R = r exp / r sim ) . 
In this way one can reduce the bias introduced by uncertainties in understanding these 
effects. 
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Super-Kamiokande, Kamiokande, 1MB and Soudan experiments 0, |J found the 



ratio R to be lower than expected, with fewer z/ M 's. Frejus and Nusex experiments ||T0| , [TT 
have not found any anomaly in R. 

If neutrinos have mass, one type of neutrino can oscillate into another. The probability 
(P) for atmospheric neutrino oscillations, in a two neutrino mixing scenario, is given by fll2"| : 

P = sin 2 29 sin 2 ( 1.27^ Am 2 ) , (1) 



E 

where 9 is the mixing angle, Am 2 is the difference of the squared mass for two types of 
neutrinos in eV 2 , L is the distance traveled by the original neutrino in km and E is the 
neutrino energy in GeV. 

The Super-Kamiokande analysis |L] constrains Am 2 from 5 x 10 -4 to 6 x 10~ 3 e^ 2 and 
9 to greater than 0.82 at 90% confidence level. The most probable solution is Am 2 = 
3.5 x 10~ 3 e\^ 2 and sin 2 29 = 1 (full mixing). Also the oscillation of a muon neutrino into a 



tau neutrino is favored over an oscillation into a sterile neutrino |T3 . 

If the atmospheric neutrino result is joined with the solar and short baseline beam mea- 
surements, they indicate evidence of a non-zero Am 2 . However, new neutrino oscillation 
experiments are needed to precisely measure Am 2 as well as to decide among experimental 
results that are in disagreement. As an example, the Super-Kamiokande [|l[] results barely 
overlap with Kamiokande results. 

As the detector area and volume covered by current or proposed high-energy astrophysical 
neutrino telescopes are large compared to underground detectors, it is possible for these 
experiments to contribute to the understanding of atmospheric neutrino oscillations. In this 
paper we investigate the possibility of AMANDA [|14j] and IceCube [^5j detectors measuring 
atmospheric neutrino oscillations. Similar analyses would apply for the ANTARES \IE and 
NESTOR M detectors. 



II. MEASUREMENT OF ATMOSPHERIC OSCILLATIONS 

Figure |] shows the neutrino flavor survival probability (Ps) versus the ratio of the distance 
traveled by the neutrino and its energy (L/E). This probability is given by 1 — P^ where 
P^ is the neutrino oscillation probability given in Equation [l|. 

As one can experimentally determine the number of events with energy E as a function 
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of the zenith angle (9z), a plot of the survival probability versus cos 9z provides a more 
convenient indication of what could be accomplished by neutrino telescopes. 

Equation [I] shows that the larger the distance the neutrino travels, the greater the pos- 
sibility of oscillation. For energies of 10 GeV and above only upward going neutrinos, those 
which come from below the detector after crossing the Earth, will have a significant proba- 
bility of oscillating. 

We define 9z as the angle from the detector axis (defined by a line from the detector to 
the center of the Earth) to the direction of the particle arrival at the detector ||33|| . For a 
detector located at a distance Rd from the center of the Earth, the distance traveled (L) by 
the neutrino will be: 

L — Rd cos 9 Z + sJrI -R 2 d + R 2 d cos 2 9 Z (2) 

where R u is the radius from the Earth's center at which the neutrino is produced, R v ~ 
R® + 15 km, where R® is the Earth radius. This equation reduces to L = Rd + R v ~ 2i? e 
for cos 9z = 1 (vertical upcoming neutrinos), to L = J ' R 2 V — R 2 , for cos 2 9z = (horizontal 
neutrinos), and L = R v — R d ~ 15 km for cos9z = —1 (vertical downgoing neutrinos). 

At the energies we are considering (10-100 GeV) only upwards going neutrinos will have a 
significant probability of oscillating. Thus the neutrino production point in the atmosphere 
will not be significant when compared to the length traveled through the Earth. 

Figure Q shows cos 9z versus the survival probability as a function of energy. We assume 
that the detector is 2 km deep in the Earth (ice or water). It can be seen that it is easier to 
detect oscillations using lower energy neutrinos. Therefore to measure neutrino oscillations 
the telescope energy threshold cannot be too high and the energy resolution must be good. 
However, for high-energy neutrino telescopes the sensing elements are placed far apart to gain 
detector volume and the trigger energy threshold is usually set high to avoid atmospheric 
muon and neutrino background. The detector design has to be optimized, if it is to be used 
for both of these two types of observations. 

Both Figure [l] and |] assume full mixing and Am 2 = 3.5 x 10~ 3 eV 2 . In Figure ||] we show 
the same dependencies as in Figure |], but for different values of Am 2 . From this figure one 
can see that, depending on the energy threshold, it may be possible to better constrain Am 2 
with a neutrino telescope. To understand whether or not this measurement can be done, 
one has to take into account the atmospheric neutrino flux, the rate of contained events and 
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FIG. 2: Neutrino survival probability versus cosine of the zenith angle for different neutrino energies 
as labeled. Full neutrino mixing and Am 2 = 3.5 x 10~ 3 eV 2 is assumed. The probability is given 
for an idealized detector 2 km deep in the Earth. It indicates that a neutrino telescope needs a 
low energy threshold (less than 30 GeV) to be able to measure atmospheric oscillations well. 

the energy resolution of the detector. 

Figure |] shows the neutrino oscillation pattern for 20 GeV neutrinos superimposed by 
other possible explanations which are compatible with the Superkamiokande |T| results. 
Measuring this oscillation pattern can allow to distinguish among different scenarios for the 
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FIG. 3: Neutrino survival probability versus cosine of the zenith angle for different neutrino energies 
(as labeled in figure Q). The probability is given for a detector 2 km deep in the Earth. Full mixing 
is assumed and Am 2 = 1.5 x 1CP 3 , Am 2 = 5.5 x 1CP 3 eV 2 as labeled. One can see that even for 
larger Am 2 the neutrino telescope energy threshold has to be equal or less than 30 GeV. 

muon neutrino deficit. 



III. ATMOSPHERIC NEUTRINO FLUX 



To determine the number of events expected to be detected, one has to convolve the 
survival probability with the atmospheric neutrino flux. We determine this flux versus cos 9z 
from the flux calculated by Volkova and compare it with other calculations [Tj^ [20f1 . 

Volkova derives the atmospheric neutrino flux from the decay of light mesons (K, tt) and 
/i's and from the decay of short-lived particles (prompt decay) which mainly consist of charm 
particles. The latter will only be significant at higher energies (around a PeV) but we include 
it for completeness. 

The differential energy spectra of atmospheric muon neutrinos from light mesons and 
muons can be approximated by MM: 
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FIG. 4: Four theoretical scenarios that are consistent with the Super-Kamionkande results. (1) 
Solid line shows best fitted standard neutrino oscillations. Neutrino survival probability versus 
cos#z for 20 GeV neutrinos crossing through the Earth (L = Earth Diameter). Full neutrino 
mixing and Am 2 = 3.5 x 10 -3 eV 2 is assumed. (2) Dashed line is effect of additional dimension 
with characteristic radius of about a micron j8). (3) Dotted line is effect of a decaying neutrino 0]. 
(4) Dashed dotted line is the effect of quantum decoherence |j] . 

= 0.48£- 4 - 04 (^(# z ) + 0.89E K (6 Z )), £ > 5.4 x 10 5 GeV, 

where E n (8z) and Ek{6z) are the ir and K critical energies, E is the neutrino energy and 
the spectra is given in units of neutrinos per cm 2 sec sr GeV. The critical energy is the one 
for which the probability for a nuclear interaction in one nuclear mean free path equals the 
decay probability in the same path. It depends on the zenith angle since the atmosphere 
density varies with depth and horizontal events cross more dense regions than the vertical 
ones. 

The spectra for prompt neutrinos can be approximated by 
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FIG. 5: Atmospheric vertical and horizontal (as labeled) muon neutrino energy spectra based on 
[18]. Also shown the contribution from K, tt and /i's and from prompt decays. 
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In Figure |5] we show the vertical and horizontal neutrino energy spectra and the contribu- 
tions from K, 7i and /x's and from prompt decays. One can see that the prompt contribution 
is negligible at lower energies. In Figure ^| we expand the lower energy region and plot the 
total flux for vertical (cos 9z = 1) and horizontal (cos^ = 0) spectra. As expected the 
number of horizontal neutrinos is slightly higher. 

We compare this flux with that obtained by |L9 1 and [ 20| . The largest discrepancy between 
these calculations is of about 15% (see Figure 14 of ref. ll| and Figure 7 of ref. p0|). In the 
energy range of interest to our work, the Volkova spectrum is an underestimate compared 
to these other spectra. We will use the atmospheric neutrino flux based on the Volkova 
spectra. This is a more conservative spectrum and having a larger flux will only enhance 
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FIG. 6: Atmospheric vertical and horizontal (as labeled) muon neutrino energy spectra. Note that 
it is plotted per km 2 instead of cm 2 as in the previous figure. The flux of horizontal neutrinos is 
higher than vertically down going neutrinos since the muons have more opportunity to decay. 

the possibility of measuring neutrino oscillations. Also important is to compare the slope 
of these spectra and Figure 7 of ref. PU| ) shows that for energies between 10 and 100 GeV 
the difference between the slope of the Volkova spectrum and of the spectrum determined 
2D | or JTP| is at maximum 1%. 



in 



We now proceed to convolve the survival probability for upward going neutrinos shown 
in Figure @ by the flux obtained from the above analysis. This is shown in Figure |7|. 
As expected, neutrino oscillations are enhanced at lower energies. 



IV. CONTAINED EVENTS 



When muon neutrinos undergo a charged current interaction, they produce a muon. This 
will propagate losing energy and eventually comes to rest and decays. If the entire event 
happens inside the detector (a "contained" event) it is much easier to reconstruct the event 
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FIG. 7: Muon neutrino survival probability weighted by atmospheric muon neutrino energy spec- 
trum. As expected, neutrino oscillations are enhanced at lower energies by a combination of higher 
flux and more rapid oscillation probability. 

than if it is not fully contained. 

The probability that a neutrino suffers a charge current interaction is given by: 



-^conv ^^cc^ (^) 

where n is the number density of nucleons of the medium transversed by the neutrino, a cc 
is the charged current cross section for a neutrino nucleon interaction and I is the distance 
traveled by the neutrino inside the detector. We determine a cc according the CTEQ4-DIS 



distributions, as described in | 21| . For energies below 10 GeV a correction on the cross 



section due to quasi-elastic and resonant effects |22| should be included. As our analysis 
starts with energies above 10 GeV we do not include these corrections. 
The number of contained events (A^nt) is therefore given by: 



N cont = P conv v3ffiA = na^^V, (6) 

where y? e is the flux of upward going neutrinos, A is the detector area and V is the detector 
volume for contained neutrino interactions. 

The neutrino flux determined in the previous section will suffer a minor attenuation when 
going through the Earth. This attenuation is not significant at these low energies but we 
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will include it for completeness. The differential flux is given by: 



dtp 

- = -na^ (7) 

where tp is the atmospherical neutrino flux, x is the distance traveled by the neutrino. The 
atmospheric neutrino flux after transversing the Earth will be given by: 

(p @ = ip oe - 1 na ~ dx ~ ^ oe - n ® <TccL ~ ^ e- 2ne<Tccfl ® cos 9z (8) 

where is the Earth radius and <^ is the initial neutrino flux, the approximations are for 
constant density and negligible detector depth, respectively. 

The volume of the detector available for the neutrino to interact and produce a contained 
muon depends on the muon range {Rp)- Making the approximation that the detector has a 
cylindrical shape, the effective volume V will be given by: 



V = hxD\ arcsin L 1 - ^| sin 2 6 Z j (l - ^| C os0 z |) (9) 

where h is the detector height and D d is the detector diameter. 

For muons with energies of tens of GeV, the average muon range can be determined 
analytically by the equation [23fl: 



fl »=H i+ H (io) 

where a and b are respectively ionization and radiation loss parameters and Eq is the initial 
muon energy. The average E is equal to the neutrino energy minus the average energy loss 
in a charged current interaction. The fractional energy loss is given in |Q and for neutrinos 
energies between 10-100 GeV is 0.48. At these energies a and b can be approximated to 
constant values where a = 2 x 10~ 3 GeVcm 2 / g and b = 4 x 10 _6 cm 2 /g [J2[j|j . 

Figure || shows the expected number of contained events per year for an idealized detector 



of the volume proposed for IceCube JT5[ with height and diameter of a kilometer. The 
neutrino energy is fixed at 20, 40 or 60 GeV and cos#z divided in bins with a width of 0.05. 
Each bin corresponds to about 0.3 sr. The same is shown for the idealized detector with the 
size of AMANDA-II 200-m diameter and 1-km height. 

As up to now we are assuming an ideal detector, these figures show that the volume of both 
AMANDA II and IceCube detectors are quite sufficient to measure atmospheric neutrino 
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FIG. 8: Annual contained events in idealized IceCube and AMANDA (as labelled) in 0.05 cos 9z 
bins versus cosine of the zenith angle. The neutrino energy is fixed as labelled. 

oscillations. The key parameter for an idealized detector is that the energy threshold has to 
be less than 30 GeV. The oscillation pattern for energies below this threshold can be seen 
when measuring the number of contained events versus the neutrino arrival direction. 

However, instrumental effects have to be taken into consideration. Before showing the sen- 
sitivity region for the Am 2 and sin 2 29 parameters that can be acheived with both AMANDA 
II and IceCube detectors, we will include these uncertainties in our analysis. 



V. INSTRUMENTAL EFFECTS AND LIMITATIONS 

In this section we consider the degradation of the measurement that results from instru- 
mental effects and limitations. Two major issues are (1) angular and (2) energy reconstruc- 
tion and resolution, since the physics manifests itself as a function of L/E and we have L 
as a function of cos 9z ■ 

There are two kinds of analysis that can be done. One measures the number of contained 
events with fixed energies (which is easy to visualize by measuring the number of events 
versus cos#z). The other measures the number of events versus the full L/E spectrum. 

To detect an oscillation pattern as a function of L/E requires collecting enough events 
and determining L/E to sufficient accuracy. From observations of the resulting muons and 
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hadronic shower one estimates the incident neutrino arrival angle 9z and energy E. The 
distance L is determined from the angle 9 (see Equation |2|). The fractional error in L/E is 
given by: 

where E v is the neutrino energy and in the approximation that L = 2i? e cos#z which is 
reasonably accurate for upward going events. For the next step we assume that the angle 
and energy correlation is small. This is likely to be only partially correct but will show 
capabilities. 



A. Angular Reconstruction and Resolution 

The muon does not travel and thus point in the same direction as the incident neutrino 
except when averaged over many events. This difference occurs because there is some ex- 
change of energy (inelasticity) and momentum to the nucleon when the neutrino converts 
into a muon. The inelasticity of the interaction and the spread in pointing are related al- 
though not as simply as one would like. Figure || shows the distribution of the angle 9 
between the muon and neutrino directions for several relevant energies. 

Characterized by one parameter, the RMS spread in direction between the incoming 



neutrino and the outgoing muon direction is given by the relation [25 



% (12) 



where the angle is given in radians. 



The RMS for 9 versus neutrino energy is shown in Figure [10]. This relation is valid for 
neutrino energies between 10 and 3000 GeV. Figure |It] shows that for neutrino energies of 
a few tens of GeV, the angular RMS varies from about 8 to 17 degrees. The RMS is spread 
enough that the direction of the neutrino will only be known by the average of the angular 
distributions of many events. 



Figure |l l (extracted from [24|) shows the mean inelasticity < y > (where 
y = (E v — E^/Ev) versus neutrino energy. In a charged current interaction < y > 
is almost constant between 10 and 100 GeV with a value of 0.48. Figure |9] shows that this 
corresponds to an average scattering angle of about 15 degrees for a 20 GeV neutrino. 
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FIG. 9: Angle distribution for several energies (as labeled). 9 corresponds to the angle between 
the neutrino and the muon directions. 

The muon direction is determined by measuring the arrival time of the Cerenkov light 
at each phototube. A muon of energy E will have a mean range of 5E m/GeV. If timing is 
done to 8t, the direction of the muon can be measured to about 

_„ V2c6t rad GeV 0.0855t rad GeV 

ou = = = 

5E m E ns 

4.86°5t GeV . . 

— ^ — 13 

E ns v ; 

For approximately upward moving muons the direction can also be determined by timing 
arrival of photons along a vertical string of tubes and comparing it to the speed of light. This 
technique measures cos Qz directly and with reasonable parameters to an accuracy better 
than 0.1. 

In the energy range that we are interested in (10 < E < 100 GeV) the angular resolution 
is therefore dominated by the neutrino-muon scattering distributions. In our analysis we 
will consider the average scattering angle. For more precision a Monte Carlo analysis can 
be done, but we consider the average scattering angle precise enough for the purpose of 
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FIG. 10: Angle distribution as a function of neutrino energy. Orms corresponds to the RMS angle 
between the neutrino and the muon directions. 

predicting the possibility of measuring neutrino oscillations. Also, Figures ^ and |3| show 
that to observe a Am 2 effect, observations do not have to be made near the horizon where 
the path length changes rapidly with angle. Thus the intrinsic spread in path length due to 
angular error is sufficiently small. 

One can therefore expect that for energies above about 10 GeV the angle of the incoming 
neutrinos can be determined sufficiently accurately for neutrino oscillation studies. 

B. Energy Reconstruction and Resolution 

We now include uncertainties on the energy estimation of the incident neutrino. An ideal 
detector, that is equally and fully sensitive to all energy deposited in the detector, could 
reconstruct the incident neutrino energy very precisely for all fully contained events. 

If all the energy from the neutrino went into the muon, then one could exploit the range 
energy relation and reconstruct the muon and thus the neutrino energy. 
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FIG. 11: Mean inelasticity versus energy. < y >= (E u — E^)/E v . Figure extracted from [24| 
We can write the formula for the error in E v as 
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where E u , E^, and .E^ are the neutrino, muon, and hadronic energies, and y = E h /E u is the 
inelasticity of the interaction. Note E v = E^ + Eh, and 1 — y = E^/ ' E v . 

The muon energy is determined by measuring its range - path length - which is on average 
5 meters per GeV. Thus a determination of the path length to an accuracy of 5 meters 
corresponds to an error of 1 GeV independent of the energy of the muon. The fluctuations 
in range have a slight energy dependence. 

One has to take the inelasticity in the conversion of neutrino to muons into account. 
Figure (extracted from ||24j| ) shows that for energies between 10 and 1000 GeV the muon 
carries about 52% of the incident neutrino energy. The inelasticity y ranges nearly uniformly 
for the energy range of interest here. At much higher energies it decreases asymptotically 
to about 20% inelasticity. 

Thus to reconstruct the neutrino energy the detector must estimate both the muon energy 
(which is given by the muon range measurement and therefore approximately independent 
of the energy of the muon) and the inelastic portion. Roughly, the energy resolution is 
composed of a component which is nearly linearly dependent on energy and a component 
which is nearly independent of the energy. The muon range is estimated by the observations 
by optical detectors of the Cerenkov radiation as it travels through the ice (or water). In 
general, the total light emitted is proportional to the energy deposited. The muon is a line 
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FIG. 12: Contained events in IceCube and AMANDA (as labelled) volume in 0.05 cos 9z bins. The 
neutrino energy is fixed as labelled and the effect of 10 and 20% energy resolution is shown. 

source of light while the other interaction products are generally more localized and near 
the point of neutrino interaction. 

The errors in calibration and fluctuations in the fraction of the energy transferred to 
the muon and to the other products of the interaction result in an error which is roughly 
proportional to the energy of the incident neutrino. This uncertainty actually decreases 
somewhat from very low energies to higher energies. 



We can therefore rewrite Equation [14| as the standard deviation of a function which 
depends on the neutrino energy: <j\ = x 2 E 2 + (SE) 2 where a; is a constant between and 
1 and 5E ~ lGeV, assuming no correlation between these terms. The energy resolution is 
therefore a fraction of the energy value plus a constant. We give our results assuming 0, 10, 
and 20% energy resolution. 

The number of contained events in IceCube and in AMANDA-II assuming the above 



energy resolution are shown in Figure [L2[ What one finds is that the standard deviation 
must be less than about 20% of the energy (or FWHM less than about 50%) in order to 
make quality oscillation observations. 

The energy threshold is also important. To be able to observe the oscillation pattern the 
energy threshold can not be much above 20 GeV. 
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FIG. 13: Flavor survival probability versus Ej cos Oz- Valid for upwards going neutrinos, where 
the distance traveled can be given in Earth radius units. The dashed line corresponds to the effect 
of 10% and the dotted dashed to 20% energy resolution. 

Figure |T3| shows the survival probability versus E/L (or E / cos 6^) including 10 and 20% 
energy resolution. In this kind of analysis, vertical upwards events (cosOz = 1) are the best 
for measuring neutrino oscillations. 

IceCube [IS| is planned to have 81 strings in a km 3 (125 m between each string). Each 



string will have 60 PMTs with 16 m between each optical module. AMANDA-II has 4 
inner strings with 20 PMTs in each of them; they are located in a cylindrical shape of 60 m 
diameter and 400 m length. There are also two outer set of PMTs; one has 6 strings with 36 
PMTs per string and 120m diameter (same height). The outermost set has 9 more strings, 
700 m long and instrumented every 14 m (50 PMTs per string). The diameter of this outer 
string array is 200 m. 

Both AMANDA and IceCube were designed to detect higher energies neutrinos. In both 
experiments, for atmospheric neutrino oscillation measurements, the strings and phototubes 
are too far apart to acheive the necessary energy resolution. Also the detector energy 
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threshold is set too high. This is true for a fixed energy analysis (which would require a 
energy threshold around 20 GeV) and for the full energy spectrum (10-100 GeV) analysis. 

For atmospheric neutrino oscillation measurements the current design proposed for Ice- 
Cub e has two major shortcomings: 

(1) Rejection of confusing events is too poor. The individual strings of optical modules 
are too far apart to guard neighboring strings against muons coming near a string and 
appearing to be a contained event or multiple muons depositing energy/light near a string 
and simulating a contained event. The rate of downward muons and muon bundles is about 
a million times the contained event rate. 

(2) The vertical optical module spacing is also too great (by a small factor) to provide 
sufficient energy resolution even for the vertical going contained events. One can estimate 
that with 16-m spacing that the range of the muon can be measured to about 8 to 10 meters 
RMS giving an error of about 2 GeV RMS for the muon energy. However about half the 
energy of the original neutrino energy is deposited in the hadron shower but in a very wide 
distribution which ranges basically from to 100%. Thus it is necessary to determine how 
much energy is deposited inelastically, which can be done by measuring how much light comes 
from the hadronic/electronic shower to the optical module and then estimating how much 
energy was deposited by determing how far it was from the shower to the optical modules. 
With a single string the vertical distance of the main energy deposition can be determined, 
though somewhat too poorly, but the distance away from the string is much more poorly 
constrained. For example, if the distance to the inelastic energy deposition is determined to 
2 to 3 meters out of a typical distance of 10 meters, then the error in estimating the distance 
translates into an error of 40% or more in reconstructing the inelastic energy. There is 
also an intrinsic spread in the amount of light received due to the statistical fluctuations in 
transit to 16-m separated optical modules. The net result is that the total energy resolution 
for vertical going contained events is too poor by a significant amount and is worse for other 
angles. 

The current version of AMANDA-II also lacks the resolution necessary for atmospheric 
neutrino oscillation measurements. At the present the trigger system and experimental 
procedures are set for measuring higher energy neutrinos with a high energy threshold (about 
50 GeV) and discriminate against contained events. 

AMANDA-II is closer to the necessary density of strings and optical modules than Ice- 
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Cube, but the vertical spacings of 20 meters, 11 meters, and 14 meters and average horizontal 
spacing of 40 meters is just too large to have the required energy resolution. The set of six 
strings with 11-meter vertical spacing would be nearly adequate but need additional infras- 
tructure in terms of other strings nearby, calibration of the optical module response and the 
appropriate triggering and data processing software. Based upon the optical properties of 
the ice and optical modules in situ and the performance of the detectors, we estimate that 
a spacing of order 5 to 10 meters would be adequate. 

Similar arguments hold when it comes to utilizing AMANDA-II and IceCube for future 
long baseline neutrino oscillation experiments [p7|| . If the beams originate from CERN 
or Fermilab, latitudes about 46° and 43° north respectively, then a typical angle of the 
incoming neutrino and produced muon is about 45° to the zenith and there is a noticeable 
spread in the muon directions. For that range of angles, the optical horizontal separation is 
as important as the vertical separation. (Also note that the rough 45° zenith angle makes 
the effective vertical separation v2 or 41% greater.) Though long-baseline experiments have 
the advantage of timing to reduce backgrounds and thus enable lower energy thresholds set 



by the experimental trigger [27], the angles for likely beams lowers the energy resolution 
substantially. The neutrino beam effective volume for IceCube would be quite low due to 
the large spacing between strings. 



VI. TOWER DETECTOR CONFIGURATION 



It is clear that measuring various values of L and E provides a very good and self- 
consistent test of neutrino oscillations. However, it is possible to test and observe neutrino 
oscillations for a fixed distance, providing sufficient range of energy and number of neutrinos 
can be observed. 

Figure [T4| shows the muon survival probability versus energy E for upwards going neu- 
trinos (L = 2R 9 ) superimposed by non-standard explanations. 

We now consider explicitly a configuration in the shape of a long tower. This tower 
is oriented vertically both to get the maximum path length L ~ 2R 9 and because that 
is operationally the most reasonable configuration for water and ice-based detectors. The 
tower geometry has its height much longer than its diameter so that its acceptance is near 
vertical (near direct upward and downward going directions) and its solid angle is quite 
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FIG. 14: Muon neutrino survival probability versus neutrino energy for upwards going neutrinos 
(L = 2R(q). Again four theoretical scenarios that might account for the observed effect are: (1) 
Standard neutrino oscillations are shown by the solid line. Full neutrino mixing and Am 2 = 
3.5 x 10~ 3 eV 2 is assumed. (2) Dashed line is effect of additional dimension with characteristic 
radius of about a micron ||. (3) Dashed Dotted line is effect of a decaying neutrino J?J]. (4) Dashed 
line is the effect of quantum gravity decoherence Q. 

limited. The limited solid angle means that the distance is effectively constant (due to the 
slow change in cos 9 near and 180 degrees). Thus the blurring due to the change in L adds 
very little to the blurring caused by the energy resolution. 

We consider the hypothetical configuration with optical detectors at 5-m to 10-m spacings 
along a string 1-km long. The close optical detector spacing improves the energy resolution 
both in terms of determining the muon range and in the energy deposited by the other 
interaction products. The tower detector consists of four such strings embedded in a larger 
detector, e.g. AMANDA-II. The larger detector acts as an after the fact additional veto 
and, when appropriate, provides additional information in constraining the event and its 
energy. 



21 




Energy (GeV) 



FIG. 15: Annual contained events for tower configuration including a 10 and 20% energy resolution. 

Figure [15| shows the number of contained events per year versus energy for vertical going 
neutrinos. In order to get the atmospheric muon neutrino flux at lower energies, we extrap- 



olate the Volkova |L8[ spectrum used in section |T| to 5 GeV. We compared this flux with 
the one calculated in [^0| and they are in good agreement. 

This figure shows the expected event rates as a function of energy with perfect energy 
resolution and including a 10 and 20% energy resolution effect. 

Previous figures in this section assume the measurement of the neutrino energy through 
measuring the muon energy and hadronic energy. Another approach to observing neutrino 
oscillations is through the measurement of the muon range, without the information on the 
hadron energy and reconstruction of the neutrino energy. Such a measurement would have 
the advantage of a clean estimate of muon energy. Since the muon carries a variable fraction 
of the energy of its parent neutrino, this is equivalent to having a random error in the 
neutrino energy estimate with a bias to the low energy side. This bias is not overwhelming 
because the atmospheric neutrino energy spectrum is steeply falling and the number of 
muons observed at energy is dominated by parent neutrino energies which are just above 
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FIG. 16: Ratio of the number of events with oscillations over the number of events without oscil- 
lations versus the neutrino energy and versus the muon energy (GeV). 

the muon energies rather than those from much higher energies. 



Figure 16 shows the ratio of the number of events with oscillations over the number of 
events without oscillations versus the neutrino energy and versus the muon energy. Mea- 
surement of neutrino oscillations using the muon energy information alone is still possible. 



This kind of analysis is proposed for ANTARES [25| and the comparison of Figure [IB] with 



Figure |2l] shows that our proposed tower configuration can measure neutrino oscillations 
as well as ANTARES, if using only the muon energy information. The key issue in that 
approach is to determine the ratio of upward going versus downward going contained muon 
events. 



VII. BACKGROUNDS AND SYSTEMATICS 



This paper is not meant to be a complete analysis of the Tower configuration for neutrino 
oscillation observations. More work is necessary. However, it is appropriate to outline and 
scope anticipated backgrounds and systematic errors to see both if there are any potential 
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show stoppers and which areas need further work. 

A. Atmospheric Muons 

The primary background is due to the very numerous downgoing muons produced by 
cosmic rays hitting the atmosphere. This is the same process that produces the atmospheric 
neutrinos. 

A veto against these events puts a more stringent requirement than that needed to ensure 
that an event is contained. A special trigger is required that looks for events and event 
topology inside the active area of the detector and vetoes against particle entering the 
detector volume from above and below or from the side. An AMANDA style low multiplicity 
majority logic trigger would be dominated by random coincidence from optical module noise. 
Deadtime due to the veto is an issue so that separate triggers would be necessary for separate 
science goals. One must also be concerned about self-vetoing due to photons leaking from 
the real signal because the absorption length in ice is fairly long. Overcoming this kind of 
self-vetoing requires careful timing which may in turn require after the fact processing of a 
much larger data set. 

The observation must have very good rejection of muons coming in at an angle and thus 
appearing to be a contained event with very little inelasticity. This is the primary motivation 
for locating the denser optical module strings inside the existing AMANDA-II array. The 
combination of the Tower closed-packed and dense strings and the information from the 19 
AMANDA-II strings provides good rejection of well off-vertical muons. 

B. Electron Neutrinos 

At low energies (order of 10 GeV), electron neutrino charged-current interactions could 
mimic low energy muon neutrino induced events. The fluxes are similar. Some are likely to 
be included in the data sample and a careful study would be needed to show that they can 
readily be rejected to the necessary level. Fortunately, one can test the results using only 
the higher energy events at some cost in sensitivity. 
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C. Neutral Current Interactions 



Although Super-Kamiokande favors muon neutrino oscillation into tau neutrino [13|, there 
is still the possibility that it oscillates into sterile neutrinos. If the latter is the case, muon 
neutrino neutral-current interactions are also a background to be considered. They will 
oscillate as the signal and due to an inelasticity comparable to the charged current one (see 
Figure |11]) will produce a jet of events from the nucleon recoil. This generates some oscillation 
in the energy spectrum of this background. The difficulty in distinguishing neutral-current 
and charge-current events is that the muon range and the photon range in ice are similar. 
So photons produced in the neutral current jet can reach the same optical modules as the 
muon Cerenkov light, and with fluctuations being what they are, it is easy to misreconstruct 
the vertex of the interaction and/or muon range. 

Vertex reconstruction is very difficult for short range (50 m) muons leaving the neutrino 
interaction vertex. Without good vertex, the energy resolution degrades. To get energy 
resolution of order 20% requires that the vertex is very well known. This requires careful 
timing and calibration of the optical modules. 

This background can be understood as long as the detection and reconstruction efficiencies 
are the same as for the charged current events. Count rates will be lower with realistic 
reconstruction efficiencies, which at these energies is limited to some extent by the muon 
neutral current interaction background. 



D. Tau-Neutrinos 



If there are neutrino oscillations, as other explanations (as shown in Figure |14]), the muon 
neutrino, most likely they oscillate into tau neutrinos fl3j . If this interpretation is correct, 



the tau neutrinos present a background which tends to wash out the oscillations. 

The outgoing tau from the tau neutrino charged current interactions decays to muons 
about 18% of the time. Thus the secondary muon will look like a muon neutrino charged 
current event and attenuate the oscillation pattern. This represents a slight decrease in 
sensitivity and is comparable in effect to an energy degradation. 

Another 11% of taus decays into a pion plus a tau neutrino, and the pion may decay to a 
muon or occasionally mimic a muon. Another 25% of the taus decays in to a pion, pi-zero, 
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and tau neutrino, and the pion decays into a muon. At 20 GeV, most pions interact before 
they decay. To do a full evaluation of these effects requires a simulation. The same may be 
required to determine what fraction of these taus will generate cascade backgrounds. 

The tau neutrino-interaction charged-current cross section has a threshold around 3 GeV, 
so any tau background is quite energy dependent and may cease to be important below 10 
GeV. This background could be quite noticeable degradation of the signal for this experiment 
and for several of the others discussed here. 



VIII. SENSITIVITY TO PARAMETERS AND COMPARISON WITH OTHER 
EXPERIMENTS 



Figure O shows the parameters sin 2 26 and Am 2 sensitivity region for the tower configu- 
ration described above. This is the region where one can demonstrate oscillation with 90% 
CL and a precision of 10% or better in both sin 2 26 and Am 2 . We assume that the detector 
energy threshold is 15 GeV and the energy resolution is 10%. The sensitivity region for the 
tower configuration is sin 2 26 > 0.40 and Am 2 > 1 x 10 -3 . 

A better precision can be achieved around the most probable values of sin 2 26 and Am 2 
(sin 2 26 = 1 and Am 2 = 3.5 x 10" 3 ). This is shown in Figures [18| and for different 



detector energy threshold and energy resolution. In this region the AMANDA II - Tower 
configuration expects to have a relative error on sin 2 26 of order 3% assuming the projected 
detector energy threshold of 15 GeV and 90% CL. The relative error on Am 2 is expected to 
be around 2% for the same energy threshold and confidence level. 

The sensitivity region and the relative error of the parameters were determined by apply- 
ing a maximum likelihood method |23] assuming a Gaussian distributed data set and using 
Fisher matrix coefficients. This analysis does not include the effect of systematic errors nor 
correlates the parameters sin 2 26 and Am 2 . However, this is not a problem when comparing 
our proposed configuration with other detectors, since none of them include systematic bias 
in their sensitivity analysis. 

A comparison with other experiments is shown in table |. MINOS |K| and MONOLITH 



3jfl cover the broadest parameter region. All experiments including our proposed tower 
configuration have roughly the same sensitivity. 

An important point is that the experiments that compose Table | represent three different 
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FIG. 17: sin 2 29 versus Am 2 sensitive region to demonstrate neutrino oscillation with the proposed 
tower configuration. The shaded region can be probed with 90% CL and with a precision of 10% 
or better in both sin 2 29 and Am 2 . A detector energy threshold of 15 GeV and a 10% energy 
resolution is assumed. 



techniques to measure neutrino oscillations. K2K and MINOS |30] use a controlled beam 



line and have two detectors, one close to the beam production and another one far away. 
Although the detectors in each experiment are different from each other, the fact that they 
can control the beam and therefore have a good energy and angular resolution is a strong 
characteristic of their design. 



MONOLITH |yj is a massive and dense magnetized tracking calorimeter. The size and 
density of the detector are important to increase the number of neutrino interactions and 
the number of contained events. MONOLITH will cover a large range in Am 2 since it is 
able to measure high momentum muons which in other experiments is not possible since the 
muons range out of the detector before losing most of their energy. 

The tower embedded in AMANDA-II that is proposed here and ANTARES are neutrino 
telescopes with large volumes in water or ice. Their advantage is their large size (detector 
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FIG. 18: 90% CL Estimated error in the mixing parameter sin 2 26 versus detector energy 
threshold. This precision can be acheived around the most probable value for neutrino oscillation 
(sin 2 29 = 1) assuming the most probable value for Am 2 (Am 2 = 3.5 x 10 3 ). 

volume). Both track and measure the energy of the muons from the Cerenkov light emitted 
in their passage through water or ice. 

The current neutrino telescopes as AMANDA, Baikal and future telescopes as IceCube, 
NESTOR and NEMO have been aimed at detecting astrophysical neutrinos and not designed 
specifically for neutrino oscillations. 

ANTARES, in contrast, has been designed to be sensitive to neutrino oscillation mea- 
surements. ANTARES [HJ has a smaller net volume than IceCube but has the potential to 
measure neutrino oscillations in the range of parameters indicated by Super-Kamionkande. 

The fact that the neutrino oscillation hypothesis will be tested with different techniques 
only enhances the importance of these experiments. All experiments will have a good sen- 
sitivity around the most probable values found by SuperKamiokande |IJ. 

Two other experimental possibilities are UNO and a long baseline experiment that 
would have a beam coming from either Fermilab or CERN CNGS to IceCube P7|. UNO 



28 




10 12.5 15 17.5 20 22.5 25 27.5 30 

Detector Energy Threshold (GeV) 



FIG. 19: 90% CL estimated error in the difference in the two neutrino mass squared Am 2 versus 
detector energy threshold for the proposed tower configuration. This precision can be achieved 
around the most probable value for neutrino oscillation (Am 2 = 3.5 x 1CP 3 ) assuming the most 
probale value for sin 2 26 (sin 2 26 = 1). 



TABLE I: Comparison among current or proposed neutrino experiments. 



Experiment 


Energy 


Am 2 (eV 2 ) 


sin 2 26 Estimated Estimated CL(%) 




Threshold 






Precision 


Precision 






(GeV) 






(Am 2 ) 


(sin 2 26) 




Proposed Tower 


16 


> 1(T 3 


> 0.4 


10% 


10% 


90 


K2K [29] 


1 


> 3 x i(r 3 


> 0.4 


30% 


? 


90 


MINOS g 


1 


> icr 3 


> 0.1 


10% 


10% 


90 


MONOLITH H 


1.5 


> 2 x i(r 4 


> 0.2 


6% 


? 


90 


ANTARES H 


5 


> 1(T 3 


> 0.6 


33% 


33% 


99 
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TABLE II: Time scale for current or proposed neutrino experiments. 
Experiment Start Time scale to achieve 

Data Taking Expected Results 
(Years) 

Proposed Tower 2003 1 
K2K M June 1999 3 



MINOS H 2003 2 



MONOLITH (31 2005-2006 4 



ANTARES m 2003 3 

would be a Cerenkov detector with 20 times the volume of Superkamiokande. It would 
identify electron and muon neutrino interactions and would have good energy and vertex 
resolution as well as good reconstruction efficiency (comparable to SuperKamiokande). UNO 
could do even longer baseline physics to look for CP violation. However this would be in 
a much longer timescale and UNO would have to adjust its goals according to new results 
obtained. 

The long baseline pTj that proposes a beamline from Fermilab or CERN to IceCube 



would have to take into account the fact that the angle which the beam would reach the 



detector is not favorable as described in section [VB| . Significant modifications would be 
necessary to the current design of IceCube to make the effort of sending a beam to IceCube 
work. 

Table |lj shows the quoted time scale for each experiment. K2K is the only one already 
in data taking mode. All others will take two more years to start collecting data. 

MINOS will be able to cover a broad sin 2 29 and Am 2 parameter space in a reasonable 
amount of time. Figure 4 of |30[ shows their precision in these parameters assuming sin 2 29 = 



0.7 and Am 2 = 5 x 10~ 3 eV 2 . In the same figure they show simulated results where, 
assuming the above values, they can distinguish among non oscillations and oscillations. 
The energy region shown in this figure ranges from zero to 20 GeV. The region which probes 
the oscillation pattern ranges from zero to less than 10 GeV. 

Figure ^ shows that the Tower configuration can probe the region above 15 GeV and 



determine the existence or not of the oscillation pattern. Figure [14| shows that the region 
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FIG. 20: L/E distribution expected from MONOLITH for Am 2 = 3 x 1(T 3 eV 2 compared to the 
best fit oscillation hypothesis (oscillating line) and to the best fit of the neutrino decay model, 
(extracted from MONOLITH proposal plj. 



which can probe the oscillation hypothesis versus non-standard scenarios that can account 
for the neutrino deficit, is above about 10 GeV. So although MINOS can probe very well the 
current parameters of neutrino oscillation, the tower configuration will be able to distinguish 
among different scenarios that can account for the neutrino deficit. Even in the case that 
the oscillation hypothesis holds, the Tower configuration and MINOS will complement each 
other results since they are probing different energy ranges in about the same timescale. 

Although MONOLITH can also probe a broad parameter space region, it can be con- 
sidered as another generation of neutrino experiments given their timescale for data taking. 



Figure ^0| (extracted from their proposal [31]) shows that they can also probe the existence 
of an oscillation pattern well above 10 GeV. In this way they will also be able to distinguish 
among different models which can account for the SuperKamiokande result. However their 
timescale is far beyond MINOS, K2K and the Tower configuration. 

K2K also has good sensitivity for the oscillation parameter space and its timescale is 
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FIG. 21: ANTARES estimated number of reconstructed events versus E/L from muon neutrinos 
undergoing oscillations with the parameters sin 2 9 = 1 (maximal coupling) and Am 2 = 0.0035 eV 2 
(points) and with no oscillations (histogram) for three years of observations. (Extracted from 
Antares proposal (J26|) . 



comparable with both MINOS and the tower configuration. 

ANTARES will also be able to test the oscillation hypothesis, since they will be able to 



go to higher energies. Figure |2l] shows the potential anticipated observations for 3 years by 
ANTARES ||28|| . Their time scale is longer than the Tower configuration since they will need 
3 years of data taking. Also the plot of their predicted results is based on the assumption 
that the atmospheric neutrino angular resolution can be measured to about 3 degrees. We 
estimate that the intrinsic spread between the observed muon and incoming neutrino is of 
order 10 degrees. 

The comparison made above shows that AMANDA-II modified to include the tower 
configuration has a strong case for probing the neutrino oscillation hypothesis. Its sensitivity 
can assure a good constraint in the current allowed region for neutrino oscillations and will 
achieve its goal in a short time scale when compared to other experiments. 
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One could argue that low energy physics should be done with lower energy detectors 
and this is not the strong suit of IceCube or AMANDA due to the photon scattering and 
low phototube density. The argument here is that with some increase in phototube density 
then the Cerenkov water detectors can be very effective for atmospheric neutrino oscillation 
measurements because they can bring large volumes to bear. 

IX. SUMMARY 

High energy neutrino telescopes may be adapted to provide very powerful observations 
of neutrino oscillations and/or some of the alternatives that might explain the up/down 
neutrino asymmetry observed by Super-Kamionkande and other anomalies. 
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Oz = degree corresponds to upward going neutrinos and 9z = 180 degrees corresponds to 
downward going neutrinos. 



Although Volkova [18| uses a lower limit of 10 2 GeV for this approximation, the extrapolation 
down to 10 GeV is in very good agreement with the flux for lower energies which can be found 
in table II of her paper. 
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